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BESCHREIBUNG 

Verfahren zur Beurteilung der Zuverlassigkeit technischer 

Systeme 

TECHNISCHES GEBIET 

Die vorliegende Erfindung bezieht sich auf das Gebiet der 
Zuverlassigkeitsanalyse technischer Anlagen. Sie geht aus 
von einem Verfahren zur quantitat iven Schatzung der Zuver- 
lassigkeit eines technischen Systems nach dem Oberbegriff 
des Anspruchs 1. 

STAND DER TECHNIK 

Ein solches Verfahren ist aus dem Artikel von T. L. Regu- 
linski und Y. P. Gupta^ "Reliability Cost Estimation: lyia- 
nagerial Perspectives", IEEE Transactions on Reliability, 
Band R-32, S. 276-281 (1983) bekannt. Dort wird ein Ver- 
fahren zur Schatzung f ehlerbedingter Kosten im Lebenszy- 
klus eines Systems dargelegt. Fur die Systemkomponenten 
werden je nach Fehlermechanismus verschiedene Ausfallwahr- 
scheinlichkeitsverteilungen postuliert und deren Parameter 
statistisch oder durch subjektive Expertenmeinung ge- 
schatzt. Eine verbesserte Expertenschat zung wird dadurch 
erzielt, dass fur die Ausfallrate ein oberer, mittlerer 
und unterer Schatzwert angegeben und zur Bestimmung einer 
Betaverteilung der Ausfallrate verwendet werden. Aus den 
solcherart modellierten Betaverteilungen wird fur das Ge- 
samtsystem eine eindeutige Kostenschat zung berechnet. Pro- 
blematisch ist jedoch, dass subjektive Schatzungen durch 
einen oder wenige Experten zu Korrelationen zwischen ver- 
schiedenen Ausfallraten fuhren und die Vorhersage fiir die 
Systemzuverlassigkeit auf unbekannte Weise systematisch 
verf alschen konnen . 
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Daruberhinaus ist es bekannt, dass fiir Fehlerarten mit 
konstanter Fehlerrate die Fehlerhauf igkeit , d. h. die Zu- 
falligkeit der Anzahl Fehler in einem Zeitintervall , durch 
eine Poissonverteilung beschrieben warden kann (s. z. B. 
das Lehrbuch von A. Birolini, "Qualitat and Zuverlassig- 
keit technischer Systeme", Springer-Verlag Berlin (1991), 
S. 366) . 

Perner sind Zuverlassigkeitsanalysen in tabellarischer 
Form unter den Begriffen FMEA-Analyse ("Failure/Fault Mo- 
des and Effects Analysis") oder FMECA ("Failure/Fault Mo- 
des, Effects and Criticality Analysis") Stand der Technik. 
Die Systemkomponenten werden mit ihren moglichen Ausfall- 
arten und den berechneten oder geschatzten Ausfallraten 
zeilenweise aufgelistet und die Storanf alligkeit des Ge— 
samtsystems beurteilt . 

DARSTELIiUNG DER ERFINDUNG 

Aufgabe der vorliegenden Erfindung ist es, ein verbesser- 
tes Verfahren zur Abschatzung der Zuverlassigkeit eines 
technischen Systeme anzugeben. Diese Aufgabe wird erfin- 
dungsgemass durch die Merkmale des Anspruchs 1 gelost , 

Die Erfindung besteht darin, dass fiir die Komponenten ei- 
nes technischen Systems erste Fehlerarten mit statistisch 
bekannten Fehlerraten und zweite Fehlerarten mit weitge- 
hend unbekannten Fehlerraten unterschieden werden, wobei 
fur letztere jeweils ein oberer und ein unterer subjekti- 
ver Schatzwert angegeben wird und mit alien oberen Schatz- 
werten eine erste Grenzwahrscheinlichkeitsverteilung und 
mit alien unteren Schatzwerten eine zweite Grenzwahr- 
scheinlichkeitsverteilung fur die Systemzuverlassigkeit 
berechnet werden. Es wird also eine ungefahre Kenntnis so- 
wohl der Fehlerraten als auch der Unsicherheit von Fehler- 
ratenschatzungen zur Angabe einer konservativen oberen und 
unteren Grenze der Systemzuverlassigkeit genutzt. 
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Mit dem erf indungsgemassen Verfahren wird der Einfluss sy- 
stematischer Fehler auf die Berechnung der Systemzuverlas- 
sigkeit weitgehend zurilckgedrangt , indem die Experten- 
schatzungen als untereinander korreliert behandelt werden. 

In einem ersten Ausf uhrungsbeispiel werden fur die zweiten 
Fehlerraten zusatzlich jeweils ein mittlerer Wert ge- 
schatzt und zur Berechnung einer mittleren Wahrscheinlich- 
keitsverteilung herangezogen . Durch Addition der einzelnen 
Ausf allwahrscheinlichkeiten bzw. durch Faltung ihrer Ver- 
teilungen kann eine grafische Darstellung der mittleren, 
maximalen und minimalen Systemzuverlassigkeit fiir einen 
Prognosezeitraum angegeben werden. Die grafische Darstel- 
lung ermoglicht eine einfache, intuitive Risikobeurteilung 
fur die Wahrscheinlichkeit eines Systemausf alls oder einer 
Kostenuberschreitung wegen Wartung, Stillstand und Repara- 
tur eines technischen Systems. 

In einem weiteren Ausf uhrungsbeispiel werden fur jede er- 
ste Fehlerart eine Poissonverteilung angenommen, fur ein 
gefordertes Konf idenzniveau eine nach oben und unten ver- 
schobene Grenz-Poissonverteilung bestimmt und durch ge- 
wichtete Summierung der drei Poissonverteilungen eine ver- 
breiterte Poissonverteilung der ersten Fehlerart bestimmt. 

Weitere Ausf uhrungen, Vorteile und Anwendungen der Erfin- 
dung ergeben sich aus den abhangigen Anspruchen sowie aus 
der nun folgenden Beschreibung anhand der Figuren. 

KURZE BESCHREIBUNG DER ZEICHNUNG 

Es zeigen beispielhaf t : 

Fig. 1 eine er f indungsgemasse Darstellung fur die prog- 
nostizierte Wahrscheinlichkeit W einer Still- 
standszeit kleiner oder gleich Tdown (in Stunden) 
wahrend eines Jahres; und 
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Fig. 2 



eine schematische Darstellung der erf indungsge- 
massen Berechnung der Systemzuverlassigkeit im 
Zusammenhang mit einer FMEA-Tabelle . 



zeichen 



In den 



Figuren sind gleiche Telle mlt gleichen Bezugs- 
versehen . 



WEGE ZUR AUSFUHRUNG DER ERFINDUNG 



Die Erflndung hat eln Verfahren zur quantltativen Schat- 
zung der Zuverlasslgkelt elnes technlschen Systems zum Ge- 
genstand. Unter System ist generell elne bellebige Maschl- 
ne Oder Anlage, z. B. elne Turbine oder eln Kraftwerk, 
Oder eln Herstellungs- oder Bearbeltungsprozess zu verste- 
hen. Dabel werden fur die Systemkomponenten Fehlerarten 
mlt zugehorlgen Fehlerraten angegeben und fiir das Gesamt- 
system und eln vorgebbares Zeltlntervall elne Wahrscheln- 
llchkeltsvertellung elner die Zuverlasslgkelt charakterl- 
slerenden Zufallsgrosse bestlmmt. Erste Fehlerraten werden 
durch statlstlsche Stichproben bestlmmt. Zwelte Fehlerra- 
ten werden durch subjektlve Expertenmelnung geschatzt, wo- 
bel jewells eln oberer und eln unterer Schatzwert angege- 
ben werden. Erf Indungsgemass werden mlt den oberen Schatz- 
werten eine erste Grenzwahrschelnlichkeltsvertellung und 
mlt den unteren Schatzwerten eine zwelte Grenzwahrscheln- 
llchkeit svertellung berechnet, Im folgenden werden bevor- 
zugte Ausf iihrungsbeispiele angegeben . 

Fur die zwelten Fehlerraten konnen zusatzllch jewells eln 
mittlerer Wert geschatzt und mlt den mlttleren Werten elne 
mlttlere Wahr sche in llchke it svertellung berechnet werden , 
Bevorzugt werden als mlttlere Werte die mutmasslich wahr- 
scheinlichsten Werte geschatzt- Das Ergebnls kann in Form 
einer grafischen Darstellung gemass Fig. 1 zusammengef asst 
und veranschaullcht werden. Dort sind als Zufallsgrosse 
elne ausf allbedlngte Stlllstandszelt Tdown wahrend elnes 
Jahres in Stunden gewahlt und drei Wahrscheinlichkeitsver- 
tellungen W (Tdown) 1/ 2, 3 dafur angegeben, dass die aus- 
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fallbedingte Stillstandszeit kleiner oder gleich Tdown ist. 
Dabei bezeichnen 2 eine erste oder optimistische und 3 ei- 
ne zweite oder pessimist ische Grenzwahrscheinlichkeitsver- 
teilung und 1 eine mittlere Wahrscheinlichkeitsverteilung . 

Fur die ersten Fehlerraten sollen Stichproben-Mittelwerte 
bestimmt und Verteilungsf unktionen angenommen werden und 
aus der Ungewissheit jedes Stichproben-Mittelwert s eine 
Breite oder eine Verbreiterung der zugehorigen Vertei- 
lungsf unktion bestimmt werden. Insbesondere werden fur je- 
de erste Fehlerart eine Poissonverteilung mit einem nomi- 
nellen Erwartungswert gleich dem Stichproben-Mittelwert 
angenommen, aus einem geforderten Konf idenzniveau ein mi- 
nimaler und ein maximaler Erwartungswert berechnet, und 
durch gewichtete Summierung der Poissonverteilungen mit 
dem nominellen, minimalen und maximalen Erwartungswert ei- 
ne verbreiterte Poissonverteilung berechnet. Erf indungsge- 
mass werden ein Konf idenzniveau 1-a und Gewichtungsf akto- 
ren l-2a, a und a fur die Poissonverteilungen mit dem no- 
minellen, minimalen und maximalen Erwartungswert gewahlt . 
Typischerweise wird a im Bereich zwischen 0,05 und 0,3 und 
insbesondere a=0,l gewahlt. 

Das erf indungsgemasse Modell einer verbreiterten Wahr- 
scheinlichkeitsverteilung dient dazu, die Ungewissheit in 
der Kenntnis des Erwartungswerts der zugrundeliegenden 
Wahrscheinlichkeitsverteilung durch eine Verbreiterung der 
Wahrscheinlichkeitsverteilung angemessen zu be rucks ichti- 
gen. Im folgenden wird ein Berechnungsbeispiel flir eine 
Poissonverteilung angegeben. Es sei aus Betriebserf ahrung 
o. a. bekannt, dass eine Systemkomponente insgesamt ngtat 
Ausfalle wahrend einer Beobachtungszeit t 

Stat hat . Dann ist 

ein nomineller Erwartungswert der zugrundegelegten Pois- 
sonverteilung durch Unom=nstat/tstat'*'T gegeben, wobei 
T=Zeitintervall flir die Riskoanalyse bzw. Betriebsdauer 
des Gesamtsystems . Man berechnet anhand eines Konfidenzni- 
veaus, z, B. l-a=0,9, eine untere und obere Grenze eines 
Konf idenzintervalls fur den Erwartungswert mit folgenden 
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Gleichungen (gemass Kosiow und Uschakow, Handbuch zur Be- 
rechnung der Zuverlassigkeit fur Ingenieure, S. 425 
(1979) ) : Uiow=niow/tstat^T, Pup=nup/ tstat*T mit 

niow^O, 5*X^ [a/2] (2nstat) und nup-=0 , 5*x^ [ l-a/2 ] (2nstat+2), wobei 
[,.,](--. ) die ( tabellierten) x^~Quantile bezeichnen. Die 
Faktoren 2 und 0,5 sind in Abweichung von einer gewohnli- 
chen durch die Poissonverteilung bedingt. Die 

Poissonverteilungen mit den Erwartungswerten jJiow/ Pnom und 
jjup werden mit den Gewichtsf aktoren 0,1; 0,8 und 0,1 ad- 
diert, so dass eine Verbreiterung resultiert. Erfindungs- 
gemass sind die Gewichtsf aktoren fur die beiden randseiti- 
gen Poissonverteilungen so gewahlt, dass die Summe der aus 
dem Konf idenzintervall nach oben und unten uberstehenden 
Anteile (jeweils ca . die Halfte der randseitigen Vertei- 
lungen) gerade gleich dem Konf idenzniveau selber sind. Da- 
durch ist die Berechnung der Poisson-Verbreiterung weitge- 
hend selbst konsistent . 

Fur ein hoheres Konf idenzniveau verbreitert sich das Kon- 
f idenzintervall , die Gewichtsf aktoren fur die Randvertei- 
lungen nehmen ab und die Verbreiterung fallt geringer aus. 
Falls das Zeitintervall T fur die Risikoanalyse deutlich 
langer als die Beobachtungszeit tstat ist, werden in der 
verbreiterten Poissonverteilung drei separate Maxima auf- 
treten. Es ist jedoch wunschenswert , eine monomodale ver- 
breiterte Poissonverteilung zu erhalten. Hierfur konnen 
die randseitigen Gewichtsf aktoren grosser oder kleiner ge- 
wahlt und/oder die Anzahl der zu gewichtenden Poissonver- 
teilungen vergrossert werden. 

Zur Verbesserung der Prognose der Systemzuverlassigkeit 
Oder der Systemrisikoanalyse werden ferner dritte Fehler- 
arten beriicksichtigt , deren Fehlerraten a priori oder mit 
grosser statistischer Sicherheit bekannt sind und fur die 
eine ideale Poissonverteilung angenommen wird. Die Vertei- 
lungsf unktionen der ersten und dritten Fehlerraten werden 
gefaltet und durch Addition der oberen, unteren und gege- 
benenfalls mittleren Schatzwerte der zweiten Fehlerraten 
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die erste und zweite Grenzwahrscheinlichkeitsverteilung 2, 
3 und gegebenenf alls eine mittlere Wahrscheinlichkeitsver- 
teilung 1 berechnet, Zusatzlich konnen andere fur die Sy- 
stemzuverlassigkeit relevante, nicht zufallige Grossen fur 
das Zeitintervall bestimmt und zu den Wahrscheinlichkeits- 
verteilungen addiert werden. 

Bei der Zuf allsgrosse, welche die Systemzuverlassigkeit 
charakterisiert , kann es sich urn eine Fehleranf alligkeit , 
eine Stillstandszeit , Kosten filr Stillstand, Reparatur und 
Wartung oder daraus abgeleitete Grossen handeln. Typi- 
scherweise sind die Fehlerraten Ausfallraten der System- 
komponenten . 

Schliesslich wird eine Anwendung des offenbarten Verfah- 
rens im Zusammenhang mit FMEA (" Failure/Fault Modes and Ef- 
fects Analysis") -Tabellen, FMECA (" Failure/Fault Modes, Ef- 
fects and Criticality Analysis ") -Tabellen oder daraus ab- 
geleiteten Tabellen erlautert. Erf indungsgemass werden in 
den Tabellen erste, zweite und dritte Fehlerraten unter- 
schieden, fur erste Fehlerraten St ichproben-Mittelwerte 
und Konf idenzintervalle angegeben, fur zweite Fehlerraten 
obere, untere und gegebenenf alls mittlere Schatzwerte an- 
gegeben und fur dritte Fehlerraten Erwartungswerte angege- 
ben - Insbesondere werden Poissonverteilungen der ersten 
und dritten Fehlerarten berechnet und miteinander gefaltet 
und werden die zweiten Fehlerraten in obere, untere und 
gegebenenf alls mittlere Fehlerraten gruppiert und jeweils 
separat addiert. 

Das Verfahren ist fur erste Fehlerarten in Fig. 2 veran- 
schaulicht. Fur jede Zeile der FMEA- oder FMECA-Tabelle 4 
mit einer ersten Fehlerart werden verbreiterte Poissonver- 
teilungen Wi berechnet 5, diese zu einer Gesamt- 
Wahrscheinlichkeitsdichteverteilung w gefaltet 6 und durch 
Integration eine Gesamt-Wahrscheinlichkeitsverteilung W 
bestimmt 7. Es kann auch die Wahrscheinlichkeitsverteilung 
1-W angegeben werden, dass die Zuf allsgrosse, z. B. die 
ausf allbedingte Stillstandszeit im vorgebbaren Zeitinter- 
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vail, grosser als Tdown ist. Fur jede Gruppe der zweiten 
Fehlerarten wird dieses Berechnungsverf ahren wiederholt, 
um eine obere, untere und gegebenenf alls mittlere Gesamt- 
Wahrscheinlichkeitsverteilung W oder 1-W fur das System- 
ausf allrisiko zu erhaiten. 



BEZUGSZEICHENIilSTE 



Tdown ausf allbedingte Stillstandszeit 

W (Tdown) Wahrscheinlichkeitsverteilung fur eine ausf allbe- 
dingte Stillstandszeit ^ Tdown 

1 mittlere Wahrscheinlichkeitsverteilung 

2 optimistische Wahrscheinlichkeitsverteilung 

3 pessimistische Wahrscheinlichkeitsverteilung 

4 FMEA-Tabelle, FMECA-Tabelle 

5 zeilenweise berechnete Wahrscheinlichkeitsdichte- 
verteilungen wi; verbreiterte Poissonverteilungen 

6 Gesamt-Wahrscheinlichkeitsdichteverteilung w; 
Faltung von ( verbeiterten) Poissonverteilungen 

7 Gesamt- Wahrscheinlichkeitsverteilung W . 
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PATENTANS PRUCHE 



Verfahren zur quant itativen Schatzung der Zuverlassig- 
keit eines technischen Systems, wobei fur die System- 
komponenten Fehlerarten mit zugehorigen Fehlerraten 
angegeben werden und fur das Gesamtsystem und fur ein 
vorgebbares Ze it interval 1 eine Wahrscheinlichkeitsver- 
teilung einer die Zuverlassigkeit charakterisierenden 
Zufallsgrosse bestimmt wird, wobei erste Fehlerraten 
durch statistische Stichproben bestimmt werden, wobei 
ferner zweite Fehlerraten durch subjektive Experten- 
meinung geschatzt werden und jeweils ein oberer und 
ein unterer Schatzwert angegeben werden, dadurch ge- 
kennzeichnet , dass mit den oberen Schatzwerten eine 
erste Grenzwahrscheinlichkeitsverteilung (2) und mit 
den unteren Schatzwerten eine zweite Grenzwahrschein- 
lichkeitsverteilung (3) berechnet werden. 

Verfahren nach Anspruch 1, dadurch gekennzeichnet , 
dass 

a) fiir die zweiten Fehlerraten zusatzlich jeweils ein 
mittlerer, insbesondere- ein wahrscheinlichster Wert 
geschatzt wird und 

b) itiit den mittleren Werten eine mittlere Wahrschein- 
lichkeitsverteilung (1) berechnet wird. 

Verfahren nach einem der Anspruche 1-2, dadurch ge- 
kennzeichnet , dass 

a) fur erste Fehlerraten Stichproben-Mittelwerte be- 
stimmt und Verteilungsf unktionen angenommen werden 
und 

b) aus der Ungewissheit jedes Stichproben-Mittelwerts 
eine Breite oder eine Verbreiterung der zugehorigen 
Verteilungsf unktion bestimmt wird. 
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Verfahren nach Anspruch 3^ dadurch gekennzeichnet , 
dass 

a) fiir jede erste Fehlerart eine Poissonverteilung mit 
einem nominellen Erwartungswert gleich dem Stich- 
proben-Mittelwert angenommen wird, 

b) aus einem geforderten Konf idenzniveau ein minimaler 
und ein maximaler Erwartungswert berechnet werden, 
und 

c) durch gewichtete Summierung der Poissonverteilungen 
mit dem nominellen, minimalen und maximalen Erwar- 
tungswert eine verbreiterte Poissonverteilung be- 
rechnet wird. 

Verfahren nach Anspruch 4, dadurch gekennzeichnet , 
dass 

a) ein Konf idenzniveau 1-a und Gewichtungsf aktoren 
l-2a, a und a fiir die Poissonverteilungen mit dem 
nominellen, minimalen und maximalen Erwartungswert 
gewahlt werden und 

b) insbesondere a=0,l gesetzt wird. 

Verfahren nach einem der Anspruche 1-5, dadurch ge- 
kennzeichnet , dass dritte Fehlerarten beriicksichtigt 
werden, deren Fehlerraten a priori oder mit grosser 
statistischer Sicherheit bekannt sind und fiir die eine 
Poissonverteilung angenommen wird. 

Verfahren nach einem der Anspriiche 1-5, dadurch ge- 
kennzeichnet , dass 

a) Verteilungsf unktionen der ersten und von dritten 
Fehlerraten gefaltet werden und durch Addition der 
oberen, unteren und gegebenenf alls mittleren 
Schatzwerte der zweiten Fehlerraten die erste und 
zweite Grenzwahrscheinlichkeitsverteilung (2, 3) 
und gegebenenf alls eine mittlere Wahrscheinlich- 
keitsverteilung (1) berechnet werden und 
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b) insbesondere andere fur die Systemzuverlassigkeit 
relevante, nicht zufallige Grossen fur das Zeitin- 
tervall bestimmt und zu den Wahrscheinlichkeitsver- 
teilungen addiert werden. 

8. Verfahren nach einem der Anspruche 1-7, dadurch ge- 
kennzeichnet , dass 

a) die die Zuverlassigkeit charakterisierende Zufalls- 
grosse eine Fehleranf alligkeit , eine Stillstands- 
zeit, Kosteri fur Stillstand, Reparatur und Wartung 
Oder daraus abgeleitete Grossen sind und 

b) insbesondere die Fehlerraten Ausfallraten sind. 

9. Anwendung des Verfahrens gemass eines der Anspruche 
1-8 im Zusammenhang mit FMEA-Tabellen, FMECA-Tabellen 
Oder daraus abgeleiteten Tabellen (4), dadurch gekenn- 
zeichnet, dass 

a) in den Tabellen (4) erste, zweite und dritte Feh- 
lerraten unterschieden werden, 

b) fiir erste Fehlerraten St ichproben-Mittelwerte und 
Konf idenzintervalle angegeben werden, 

c) fiir zweite Fehlerraten obere, untere und gegebenen- 
falls mittlere Schatzwerte angegeben werden und 

d) fiir dritte Fehlerraten Erwartungswerte angegeben 
werden . 

10. Anwendung des Verfahrens nach Anspruch 9, dadurch ge- 
kennzeichnet , dass 

a) Poissonverteilungen der ersten und dritten Fehler- 
arten berechnet und miteinander gefaltet werden und 

b) die zweiten Fehlerraten in obere, untere und gege- 
benenfalls mittlere Fehlerraten gruppiert und je- 
weils separat addiert werden. 
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ZUS AMMENFAS SUNG 



Die vorliegende Erfindung betrifft ein Verfahren zur quan- 
titativen Schatzung der Zuverlassigkeit eines technischen 
Systems, das besonders flir komplexe Systems mit einer 
Vielzahl von Komponenten von Nutzen ist. Fur die System- 
komponenten warden erste, zweite und dritte Ausfallraten 
unterschieden . Erf indungsgemass werden fur die weitgehend 
unbekannten zweiten Ausfallraten jeweils ein oberer, unte- 
rer und gegebenenf alls mittlerer Wert durch subjektive Ex- 
pertenmeinung geschatzt und alle oberen, unteren und gege- 
benenfalls mittleren Schatzwert getrennt zur Berechnung 
einer optimist ischen und pessimistischen Grenzwahrschein- 
lichkeitsverteilung (2, 3) und gegebenenf alls einer mitt- 
leren Wahrscheinlichkeitsverteilung (1) der Systemzuver- 
lassigkeit verwendet. Dadurch werden systematische Korre- 
lationen zwischen Expertenschat zungen berucksichtigt . Fur 
die ersten Ausfallraten werden aus Betriebserf ahrung ge- 
wonnene Mittelwerte mit einem Konf idenzintervall bestimmt 
und verbreiterte Poissonverteilungen berechnet . Zur Ge- 
samtwahrscheinlichkeit der Systemzuverlassigkeit konnen 
ferner die Poissonverteilungen von dritten Fehlerarten ad- 
diert werden, deren Ausfallraten a priori oder mit grosser 
statistischer Sicherheit bekannt sind. Das Verfahren ist 
besonders fur die Verwendung im Zusammenhang mit 
FMEA ( "Failure/Fault Modes and Effects Analysis ") -Tabellen 
(4) und FMECA ( "Failure/Fault Modes, Effects and Criticali- 
ty Analysis" ) -Tabellen (4) geeignet. 



(Fig. 1) 
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Method for assessing the reliability of technical 

systems 

The present invention relates to the field of 
reliability analysis of technical systems. It proceeds 
from a method for quantitative estimation of the 
reliability of a technical system according to the 
preamble of claim 1. 

Such a method is known from the article by T. L. 
Regulinski and Y. P. Gupta, "Reliability Cost 
Estimation: Managerial Perspectives", IEEE Transactions 
on Reliability, Volume R-32, pages 276-281 (1983). A 
method is set forth there for estimating fault-induced 
costs in the life cycle of a system. Various failure 
probability distributions are postulated for the system 
components, depending on the fault mechanism, and their 
parameters are estimated statistically or by subjective 
expert opinion." An improved expert estimate is achieved 
by specifying an upper, mean and lower estimate for the 
failure rate and using them to determine a beta 
distribution of the failure rate. A unique cost 
estimate is calculated for the overall system from the 
beta distributions modeled in such a way. However, 
there is a problem that subjective estimates by one or 
a few experts lead to corrections between different 
failure rates, and the systematic falsification in an 
unknown way of the prediction of the system 
reliability. 

Moreover, it is known that for types of fault with a 
constant fault rate the fault frequency, that is to say 
the randomness of the number of faults in a time 
interval, can be described by a Poisson distribution 
(see, for example, the textbook by A. Birolini, 
"Qualitat und Zuverlassigkeit technischer Systeme" 
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[Quality and Reliability of Technical Systems] , 
Springer-Verlag Berlin (1991), page 366). 

Furthermore, reliability analyses in tabular form under 
the terms of FMEA analysis ("Failure/Fault Modes and 
Effects Analysis") or FMECA ("Failure/Fault Modes, 
Effects and Criticality Analysis") constitute prior 
art. The system components are listed in rows with 
their possible types of failure and the calculated or 
estimated failure rates and the susceptability to 
disruption of the overall system are assessed. 

The object of the present invention is to specify an 
improved method for estimating the reliability of a 
technical system. This object is achieved according to 
the invention by means of the features of claim 1. 

The invention consists in that first types of fault 
with statistically known fault rates and second types 
of fault with largely unknown fault rates are 
distinguished for the components of a technical system, 
an upper and a lower subjective estimate being 
specified respectively for the latter fault rates, and 
a first limiting probability distribution being 
calculated for the system reliability with all the 
upper estimates and a second limiting probability 
distribution being calculated for the system 
reliability with all the lower estimates. Thus, 
approximate knowledge both of the fault rates and of 
the uncertainty of fault rate estimates is used to 
specify a conservative upper and lower limit of the 
system reliability. 

The method according to the invention largely cuts back 
the influence of systematic faults on the calculation 
of the system 'reliability by treating the expert 
estimates as being intercorrelated . 
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In a first exemplary embodiment, a mean value is 
respectively estimated in addition for the second fault 
rates and used to calculate a mean probability 
distribution. By adding the individual failure 
probabilities and by convoluting their distributions, 
it is possible to specify a graphical representation of 
the mean, maximum and minimum system reliability for a 
prediction period. The graphical representation permits 
a simple, intuitive risk assessment for the probability 
of a system failure or a cost overrun owing to 
maintenance, standstill and repair of a technical 
system . 



In a further exemplary embodiment, a Poisson 
distribution is assumed for each first type of fault, 
an upwardly and downwardly displaced limiting Poisson 
distribution is determined for a required confidence 
level, and a widened Poisson distribution of the first 
type of fault is determined by weighted summing of the 
three Poisson distributions. 



Further embodiments, advantages and applications of the 
invention follow from the dependent claims and from the 
description now following with the aid of the figures. 

By way of example: 

Figure 1 shows a representation according to the 
invention for the predicted probability W of 
a downtime of less than or equal to Tdown{in 
hours) in the course of one year; and 

Figure 2 shows a schematic representation of the 
calculation according to the invention of the 
system reliability in conjunction with an 
FMEA table. 



Identical parts -are provided with identical reference 
symbols in the figures. 



- 4 - 



99/045 



The subject matter of the invention is a method for 
quantitative estimation of the reliability of a 
technical system. A system is understood to be in 
general terms an arbitrary machine or system, for 
example a turbine or a power plant, or a process for 
production or treatment. In this case, types of fault 
with associated fault rates are specified for the 
system components, and a probability distribution of a 
random variable characterizing the reliability is 
determined for the overall system and a prescribable 
time interval. First fault rates are determined by 
statistical random samples. Second fault rates are 
estimated by subjective expert opinion, an upper and a 
lower estimate being specified in each case. According 
to the invention, a first limiting probability 
distribution is calculated with the upper estimates, 
and a second limiting probability distribution is 
calculated with the lower estimates. Preferred 
exemplary embod^iments are specified below. 

For the second fault rates, it is additionally possible 
to estimate a mean value in each case and to calculate 
a mean probability distribution with the mean values. 
The values presumed to be the most probable are 
preferably estimated as mean values. The result can be 
sximmarized and illustrated in the form of a graphical 
representation in accordance with Figure 1. In the 
figure, a failure- induced downtime Tdown during a year 
in hours is selected as random variable, and three 
probability distributions W (Tdown) 1/ 2, 3 are specified 
for the failure-induced downtime being less than or 
equal to Tdown- In this case, 2 denotes a first or 
optimistic, and 3 a second or pessimistic, limiting 
probability distribution, and 1 denotes a mean 
probability distribution. 

The aim for the first fault rates is to determine 
random sample mean values and assume distribution 
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functions, and to determine from the uncertainty of 
each random sample mean value a width or a widening of 
the associated distribution function. In particular, 
for each first type of fault, a Poisson distribution 
with a nominal expectation value equal to the random 
sample mean value is assiimed, a minimum and a maximum 
expectation value are calculated from a required 
confidence level, and a widened Poisson distribution is 
calculated by weighted summing of the Poisson 
distributions with the nominal, minimum and maximum 
expectation values. According to the invention, a 
confidence level 1-a and weighting factors l-2a, a and 
a are selected for the Poisson distributions with the 
nominal, minimum and maximum expectation values, a is 
typically selected in the range between 0.05 and 0.3, 
and in particular a=0 . 1 . 

The model according to the invention of a widened 
probability dis^tribution serves the purpose of taking 
appropriate account of the uncertainty in the knowledge 
of the expectation value of the fundamental probability 
distribution by widening the probability distribution. 
A computational example is specified below for a 
Poisson distribution. It may be assumed as known from 
operational experience or the like that a system 
component has a total of nstat failures during one 
observation period tstat- A nominal expectation value of 
the fundamental Poisson distribution is then given by 
Mnoin=nstat/ tstaii*T, . where T=time interval for the risk 
analysis or operating period of the overall system. 
Using a confidence level of, for example, l-a=0.9, a 
lower and upper limit of a confidence interval is 
calculated for the expectation value using the 
following equations (after Koslow and Uschakow, 
Handbuch zur Berechnung der Zuverlassigkeit fur 
Ingenieure, [Engineer's manual for calculating 
reliability] , page 426 (1979) ) : |iiow=niow/ tstat*T , 
|iaip=nup/ tstat *T, where niow=0 . 5*%^ [a/2 ] (2nstat) and 
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nup=0.5*x^[l-a/2] (2nstat+2) , where (...) denote the 

(tabulated) quantiles. In a departure from a 

customary distribution, the factors 2 and 0.5 are 

determined by the Poisson distribution. The Poisson 
distributions with the expectation values M-iow/IXnom and 
|Xup are added to the weighting factors 0.1; 0.8 and 0.1, 
thus resulting in widening. According to the invention, 
the weighting factors for the two edge Poisson 
distributions are selected such that the sum of the 
components projecting upwards and downwards from the 
confidence interval (in each case, approximately half 
of the edge distributions) are exactly equal to the 
confidence level itself. The calculation of the Poisson 
widening is thereby largely self -consistent . 

For a higher confidence level, the confidence interval 
is widened, the weighting factors for the edge 
distributions decrease, and the widening is 
disminished. If the time interval T for the risk 
analysis is substantially longer than the observation 
time tstat/ three separate maxima occur in the widened 
Poisson distribution. However, it is desirable to 
obtain a monomodal widened Poisson distribution. The 
edge weighting factors can be selected to be larger or 
smaller for this purpose, and/or the number of the 
Poisson distributions to be weighted can be enlarged. 

Furthermore, in order to improve the prediction of the 
system reliability or the system risk analysis, account 
is taken of third types of fault whose fault rates are 
known a priori or with high statistical reliability, 
and for which an ideal Poisson distribution is assumed. 
The distribution functions of the first and third fault 
rates are convoluted, and the first and second limiting 
probability distributions 2, 3 and, if appropriate, a 
mean probability distribution 1 are calculated by 
adding the upper, lower and, if appropriate, mean 
estimates of the second fault rates. In addition, other 
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non-random variables relevant to the system reliability 
can be determined for the time interval and added to 
the probability distributions. 

The random variable which characterizes the system 
reliability can be a susceptibility to faults, 
downtime, costs of standstill, repair and maintenance, 
or variables derived therefrom. The fault rates are 
typically failure rates of the system components. 

Finally, there is an explanation of an application of 
the disclosed method in conjunction with FMEA 
("Failure/Fault Modes and Effects Analysis") tables, 
FMECA ("Failure/Fault Modes, Effects and Criticality 
Analysis") tables, or tables derived therefrom. 
According to the invention, first, second and third 
fault rates are distinguished in the tables, random 
sample mean values and confidence intervals are 
specified for first fault rates, upper, lower and, if 
appropriate, mean estimates are specified for second 
fault rates, and expectation values are specified for 
third fault rates. In particular, Poisson distributions 
of the first and third types of fault are calculated 
and convoluted with one another, and the second fault 
rates are grouped into upper, lower and, if 
appropriate, mean fault rates and added separately in 
each case. 

The method is illustrated for the first types of fault 
in Figure 2 . Widened Poisson distributions wi are 
calculated 5 for each row of the FMEA or FMECA table 4 
with a first type of fault and convoluted 6 to form an 
overall probability density distribution W, and an 
overall probability distribution W is determined 7 by 
integration. It is also possible to specify the 
probability distribution 1-W such that the random 
variable, for example the failure-induced downtime in 
the prescribable time interval, is greater than Tdown- 
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This computing method is repeated for each group of 
second types of fault in order to obtain an upper, 
lower and, if appropriate, mean overall probability 
distribution W or 1-W for the system failure risk. 
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LIST OF DESIGNATIONS 



Tdown Failure-induced downtime 

W(Tdown) Probability distribution for a failure-induced 
downtime < Tdown 

1 Mean probability distribution 

2 Optimistic probability distribution 

3 Pessimistic probability distribution 

4 FMEA table, FMECA table 

5 Probability density distributions Wi calculated 
by rows; widened Poisson distributions 

5 Overall probability density distribution w; 

convolution of (widened) Poisson distributions 
7 Overall probability distribution W 
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PATENT CLAIMS 

A method for quantitative estimation of the 
reliability of a technical system, types of fault 
with associated fault rates being specified for 
the system components, and a probability 
distribution of a random variable characterizing 
the reliability being determined for the overall 
system and for a prescribable time interval, first 
fault rates being determined by statistical random 
samples, and second fault rates furthermore being 
estimated by subjective expert opinion and an 
upper and a lower estimate being specified in each 
case, characterized in that a first limiting 
probability distribution (2) is calculated with 
the upper estimates, and a second limiting 
probability distribution (3) is calculated with 
the lower estimates. 

The method as claimed in claim 1, characterized in 
that 

a) a mean, in particular a most probable, value 
is additionally estimated in each case for 
the second fault rates, and 

b) a mean probability distribution (1) is 
calculated with the mean values . 

The method as claimed in one of claims 1-2, 
characterized in that 

a) random sample mean values are determined for 
first fault rates and distribution functions 
are assumed, and 

b) a width or a widening of the associated 
distribution function is determined from the 
uncertainty of each random sample mean value. 

The method as claimed in claim 3, characterized in 
that 
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a) for each first type of fault, a Poisson 
distribution with a nominal expectation value 
equal to the random sample mean value is 
assumed, 

b) a minimum and a maximum expectation value are 
calculated from a required confidence level, 
and 

c) a widened Poisson distribution is calculated 
by weighted siomming of the Poisson 
distributions with the nominal, minimum and 
maximum expectation values. 

The method as claimed in claim 4, characterized in 
that 

a) a confidence level 1-a and weighting factors 
l-2a, a and a are selected for the Poisson 
distributions with the nominal, minimum and 
maximum expectation values, and 

b) in particular, a=0.1 is set. 

The method as claimed in one of claims 1-5, 
characterized in that account is taken of third 
types of fault whose fault rates are known a 
priori or with high statistical reliability, and 
for which a Poisson distribution is assumed. 

The method as claimed in one of claims 1-6, 
characterized in that 

a) distribution functions of the first and third 
fault rates are convoluted, and the first and 
second limiting probability distributions (2, 
3) and, if appropriate, a mean probability 
distribution (1) are calculated by adding the 
upper, lower and, if appropriate, mean 
estimates of the second fault rates, and 

b) in particular, other non-random variables 
relevant to the system reliability are 
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determined for the time interval and added to 
the probability distributions . 



The method as claimed in one of claims 1-7, 
characterized in that 

a) the random variable characterizing the 
reliability is a susceptability to faults, a 
downtime, costs of standstill, repair and 
maintenance, or variables derived therefrom, 
and 

b) in particular, the fault rates are failure 
rates . 



The application of the method in accordance with 
one of claims 1-8 in conjunction with FMEA tables, 
FMECA tables or tables (4) derived therefrom, 
characterized in that 

a) first, second and third fault rates are 
distinguished in the tables (4) , 

b) random sample mean values and confidence 
intervals are specified for first fault 
rates , 

c) upper, lower and, if appropriate, mean 
estimates are specified for second fault 
rates , and 

d) expectation values are specified for third 
fault rates. 

The application of the method as claimed in claim 
9, characterized in that 

a) Poisson distributions of the first and third 
types of fault are calculated and convoluted 
with one another, and 

b) the second fault rates are grouped into 
upper, lower and, if appropriate, mean fault 
rates and added separately in each case. 
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The present invention relates to a method for 
quantitative estimation of the reliability of a 
technical system, which is useful, in particular, for 
complex systems with a multiplicity of components. 
First, second and third failure rates are distinguished 
for the system components. In accordance with the 
invention, an upper, lower and, if appropriate, mean 
value are in each case estimated for the largely 
unknown second failure rates by subjective expert 
opinion, and all upper, lower and, if appropriate, mean 
estimates are used separately to calculate an 
optimistic and pessimistic limiting probability 
distribution (2, 3) and, if appropriate, a mean 
probability distribution (1) of the system reliability. 
Systematic correlations between expert estimates are 
thereby taken into account. For the first failure 
rates, mean values obtained from operational experience 
are determined with a confidence interval, and widened 
Poisson distributions are calculated. Furthermore, it 
is possible to add to the overall probability of the 
system reliability the Poisson distributions of third 
types of fault whose failure rates are known a priori 
or with a high statistical reliability. The method is 
suitable, in particular, for use in connection with 
FMEA ("Failure/Fault Modes and Effects Analysis") 
tables (4) and FMECA ("Failure/Fault Modes, Effects and 
Criticality Analysis") tables (4). 



(Figure 1) 



